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Ultracold atoms in optical lattices are a versatile tool to investigate fundamental properties of 
quantum many body systems. In particular, the high degree of control of experimental parame- 
ters has allowed the study of many interesting phenomena such as quantum phase transitions and 
quantum spin dynamics. Here we demonstrate how such control can be extended down to the most 
fundamental level of a single spin at a specific site of an optical lattice. Using a tightly focussed 
laser beam together with a microwave field, we were able to flip the spin of individual atoms in 
a Mott insulator with sub-diffraction-limited resolution, well below the lattice spacing. The Mott 
insulator provided us with a large two-dimensional array of perfectly arranged atoms, in which we 
created arbitrary spin patterns by sequentially addressing selected lattice sites after freezing out the 
atom distribution. We directly monitored the tunnelling quantum dynamics of single atoms in the 
lattice prepared along a single line and observed that our addressing scheme leaves the atoms in 
the motional ground state. Our results open the path to a wide range of novel applications from 
quantum dynamics of spin impurities, entropy transport, implementation of novel cooling schemes, 
and engineering of quantum many-body phases to quantum information processing. 



The ability to observe and control the position of sin- 
gle atoms on a surface of a solid via scanning tunnelling 
and atomic force microscopy has revolutionised the field 
of condensed matter physics [1, 2 . In few-atom systems, 
coherent control of single particles in e.g. an ion chain 
has proven crucial for the implementation of high-fidelity 
quantum gates and the readout of individual qubits in 
quantum information processing [3 . Bringing such lev- 
els of control to the regime of large scale many-body sys- 
tems has been a longstanding goal in quantum physics. 
In the context of ultracold atoms in optical lattices, a 
major challenge has been to combine degenerate atomic 
samples with single-site addressing resolution and single- 
atom sensitivity. This full control is essential for many 
applications in condensed matter physics, such as the 
study of spin impurities [4] and quantum spin dynamics 
[5l|6] within quantum magnetism, entropy transport, the 
implementation of novel cooling schemes [Zl |8] or digital 
quantum simulations based on Rydberg atoms [9 . For 
scalable quantum information processing, a Mott insula- 
tor with unity filling provides a natural quantum register 
with several hundreds of qubits. In order to exploit the 
full potential of such a large scale system for quantum 
computation, coherent manipulation of individual spins 
is indispensable, both within a circuit-based [10^ or a 
one-way quantum computer architecture [TTJ [12] . 

The quest to address atoms on single sites of an opti- 
cal lattice has a long history [7^ .13-22^ . In one dimen- 
sion, single-site addressing was accomplished optically 
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in a long- wavelength lattice [23^, in which however tun- 
nelling was completely suppressed, and using magnetic 
resonance techniques in a sparsely filled short-wavelength 
lattice 24l[25]. In two dimensions, an electron beam was 
used to depopulate sites of a Bose-Einstein condensate 
loaded into an array of potential tubes, each containing 
up to 80 atoms [26]. In this case, coherent spin manip- 
ulation was not possible and the readout was done by 
averaging over more than 100 single images. None of the 
experiments so far has shown single-atom spin control 
in strongly correlated systems together with high fidelity 
single-atom detection. 

Here we report on the achievement of this goal, by 
demonstrating single-site-resolved addressing and control 
of the spin states of individual atoms in a Mott insula- 
tor in an optical lattice. The Mott insulator provided 
us with an almost perfect initial two-dimensional array 
of atoms in the same spin state. Apart from few ther- 
mal defects, each lattice site contained a single atom in 
its motional ground state [27l [28]. Using a tightly fo- 
cused laser beam, we introduced a controlled differential 
energy shift between two atomic spin states at a given lat- 
tice site. Microwave radiation resonant with this shifted 
transition then allowed us to selectively address the spin 
of a single atom [3 |T7] with high fidelity. We thereby ob- 
tained sub-diffraction-limited spatial resolution well be- 
low the lattice spacing. By moving the addressing laser 
beam to different lattice sites and by inducing spin-flips 
in the Mott insulator, we were able to deter ministically 
create arbitrary two-dimensional spin patterns of indi- 
vidual atoms, thereby realising a scalable single-atom 
source [29^-^3^. Furthermore, we investigated how much 
our single-spin manipulation affects the motional state of 
the atoms by directly monitoring the tunnelling dynam- 
ics of single atoms after addressing them. Averaging over 
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FIG. 1: Addressing scheme, a, Atoms in a Mott insula- 
tor with unity filling arranged on a square lattice with pe- 
riod aiat = 532 nm were addressed using an off-resonant laser 
beam. The beam was focussed onto individual lattice sites by 
a high- aperture microscope objective (not shown) and could 
be moved in the xy plane with an accuracy of better than 
O.laiat- b, Energy diagram of atoms in the lattice for the 
two hyperfine states |0) and The cr~-polarized addressing 
beam locally induces a light shift Als of state |1), bringing 
it into resonance with a microwave field. A Landau- Zener 
sweep (central frequency cjmw, sweep width ctmw) transfers 
the addressed atoms from |0) to 

several snapshots after different tunnelling times, we fully 
reconstructed the characteristic spatial probability distri- 
bution of the single-atom wave function and its coherent 
evolution over more than 20 lattice sites. We were able 
to discriminate the dynamics of the atoms in the zeroth 
and in the first band and found that most of the atoms 
remained in the motional ground state after addressing. 



Experimental setup 

In our experiments, we prepared a two-dimensional 
sample of ultracold ^^Rb atoms in an optical lattice, con- 
fined in a single antinode of a vertical standing wave along 



the z direction. Two pairs of counterpropagating laser 
beams (wavelength A = 1064 nm) along the horizontal x 
and y directions provided a square lattice with period of 
<^iat = A/2 = 532 nm (for details see Ref. [28 ). Starting 
from a Bose-Einstein condensate we raised the poten- 
tial in the x and y lattice axes within 75 ms to values of 
Vx.y — 23(2) (the number in parenthesis denotes the 
Icr uncertainty of the last digit), where Ej. = /{2m\^) 
is the recoil energy, and m denotes the atomic mass of 
^^Rb. In this way, the interaction energy of the particles 
with respect to their kinetic energy was increased such 
that the system undergoes a quantum phase transition 
to a Mott insulating state |32H33] . Due to the external 
harmonic confinement, the Mott insulator exhibits a shell 
structure with fixed integer atom numbers increasing in 
a step-like manner from the outer regions of the system 
to the inner core [35] [36]. The initial state for all ex- 
periments presented in this paper was a single shell with 
only one atom per lattice site, which in our system was 
realised for atom numbers smaller than ~ 400. 

We detected the atoms using fluorescence imaging via 
a high-resolution microscope objective with numerical 
aperture of NA = 0.68. An optical molasses induced fluo- 
rescence light and simultaneously laser-cooled the atoms. 
Light-assisted collisions lead to rapid losses of atom pairs, 
such that we only detected the atom number modulo two 
on each lattice site [27l[28]. With about 5,000 collected 
photons per atom, we identified individual atoms in the 
lattice with an excellent signal-to-noise ratio. Even in 
the regions of high atomic density, we determined the 
presence or absence of an atom for each lattice site with 
> 99.5% fidelity using a special reconstruction algorithm 



Addressing single lattice sites 

In order to address the atoms in the lattice, we used 
an off-resonant laser beam focused by the high-resolution 
imaging system onto individual lattice sites (Fig.[T]). The 
laser beam causes a differential light shift of the two 
relevant hyperfine levels |0) = |F = l,mi? = — 1) and 
|1) = \F = 2, mi? = — 2) and tunes the addressed atom 
into resonance with an external microwave field at ~ 6.8 
GHz. The <j~ -polarized addressing beam had a wave- 
length of 787.55 nm, between the Di and D2 lines, in 
order to obtain a large differential light shift between 
the two hyperfine levels. For perfect <j~ -polarization, 
this 'magic' wavelength generates a light shift only for 
state and leaves state |0) unaffected. The beam 
had a diameter of 600 nm full-width at half-maximum 
(FWHM) and could be moved in the object plane over 
the entire field of view by changing its angle of incidence 
into the microscope objective with a two-axis piezo mir- 
ror. We were able to position the beam with an accu- 
racy better than 0.1 aiat using an independent calibration 
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FIG. 2: Single-site addressing, a, Experimentally obtained fluorescence image of a Mott insulator with unity filling in which 
the spin of selected atoms was flipped from |0) to |1) using our single-site addressing scheme. Atoms in state |1) were removed 
by a resonant laser pulse before detection. The lower part shows the reconstructed atom number distribution on the lattice. 
Each circle indicates a single atom, the points mark the lattice sites, b-f, Same as a, but a global microwave sweep exchanged 
the population in |0) and |1), such that only the addressed atoms were observed. The line in b shows 14 atoms on neighbouring 
sites, the images c-f contain 29, 35, 18 and 23 atoms, respectively. The single isolated atoms in b,e,f were placed intentionally 
to allow for the correct determination of the lattice phase for the feedback on the addressing beam position. 



measurement of its position together with a feedback that 
tracks the slowly varying lattice phases (see Appendix). 

If the addressing laser beam is perfectly centered onto 
one lattice site (see Fig.[l]3), the differential light shift 
is ALs/(27r) —70kHz, whereas a neighbouring atom 
only experiences 10% of the peak intensity. The result- 
ing difference in light shifts can be well resolved spec- 
trally by our microwave pulses. In order to flip the spin, 
we performed Landau- Zener sweeps (see Appendix) of 
c^Mw/(27r) = 60 kHz width and 20 ms duration yielding a 
near flat-top frequency spectrum with a maximum pop- 
ulation transfer efficiency of ^ 95%. 

As a first experiment, we sequentially flipped the spin 
of the atoms at selected lattice sites in our Mott insulator 
with unity filling and spin state |0) (Fig. [2|l). The lattice 
depths were first changed to Vx = 56Er^ Vy = 90 Er and 
Vz = 70 Er in order to completely suppress tunnelling 
even when the addressing beam locally perturbs the lat- 
tice potential. For each lattice site, we then switched on 
the addressing laser beam with an s-shaped ramp within 
2.5 ms, which is adiabatic with respect to the on-site 
oscillation frequency of ~ 30 kHz. Subsequently, a mi- 
crowave pulse with the parameters described above pro- 
duced spin-flips from |0) to |1). The addressing laser was 
switched off again within 2.5 ms, before its position was 
changed in 5 ms to address the next lattice site. For the 
image of Fig. [2^, this procedure was repeated 16 times in 



order to flip the spins along a line. Finally, a 5 ms 'push- 
out' laser pulse, resonant with the F = 2 to = 3 tran- 
sition, removed all addressed atoms in state |1). In order 
to reveal only the spin-flipped atoms, the spin states of all 
atoms were flipped by a global microwave sweep before 
the push-out laser was applied (Fig.[2]3,d). However, due 
to the finite transfer efficiency of the global sweep, some 
atoms remaining in state |0) were clearly visible in addi- 
tion to the addressed ones. To avoid this problem when 
detecting the addressed atoms, we initially transferred 
the whole sample to state |1) by a microwave sweep and 
then shone in a repumping laser that completely depop- 
ulated state |0). Then, we used our addressing scheme 
to transfer selected atoms back to |0) and subsequently 
pushed out the atoms in |1), yielding typical images as 
shown in Fig.[2]3,e,f. 

Spin-flip fidelity 

We quantified the success rate of our addressing scheme 
by again producing a series of spin-flips along the y lattice 
axis in a Mott insulator with unity filling (see Fig.[3|. The 
experimental procedure was the same as described above 
for the realisation of Fig. [2^, in which the addressed sites 
were detected as empty sites. From the reconstructed 
atom number distribution (Fig.|3|i), we determined the 
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FIG. 3: Addressing fidelity. The spin- flip probability was 
measured by sequentially addressing a series of 16 neighbour- 
ing sites along the y lattice axis (red circles in a) in a Mott 
insulator with unity filling. The red data points in c show 
the resulting hole probability po{Sx) as a function of the 
pointing offset 6x, as defined in b. Each point was obtained 
by averaging over 4 — 7 pictures (total 50 — 100 addressed 
lattice sites), taking only those sites into account which lie 
well within a Mott shell with unity filling. The displayed 
error bars show the la statistical uncertainty, given by the 
Clopper- Pearson confidence limits. The data was fitted by 
a flat-top model function (see Appendix) and yields a full- 
width at half-maximum aa = 330(20) nm, an edge sharpness 
of as = 50(10) nm, and a peak fidelity of 95(2)%. The offset 
was fixed at the 6(2)% probability of thermally activated holes 
as deduced from the neighbouring and next neighbouring sites 
(blue shaded regions in a,b and blue points in c). 



probability po{Sx) of finding an empty site as a function 
of the pointing offset 5x between the addressing beam 
and the center of the lattice site (see Fig.[3]3). We also 
investigated the effect of the addressing on neighbouring 
atoms, which ideally should remain unaffected. For this 
purpose, we monitored the probability of finding a hole 
at the sites next to the addressed ones (dark blue regions 
in Fig.|3^,b and points in Fig.js]^). In order to distinguish 
accidentally flipped neighbouring atoms from holes that 
originate from thermal excitations of the initial Mott in- 
sulator [28] , we also monitored the probability of finding 
a hole at the second next neighbours (light blue regions 
and points in Fig.|3|. As both yielded the same hole 
probability of 6(2)%, we attribute all holes to thermal 
excitations and conclude that the probability of address- 



ing a neighbouring atom is indiscernibly small. We fitted 
the hole probability po{Sx) of the addressed site with a 
flat-top model function (see Appendix), keeping the off- 
set fixed at the thermal contribution of 6%. From the fit, 
we derived a spin- flip fidelity of 95(2)%, a full- width at 
half-maximum of a a = 330(10) nm and an edge sharpness 
of (Is = 50(10) nm (Fig.jsj^). These values correspond to 
60% and 10% of the addressing beam diameter, demon- 
strating that our method reaches sub-diffraction-limited 
resolution, well below the lattice spacing. 

The observed maximum spin-flip fidelity is currently 
limited by the population transfer efficiency of our mi- 
crowave sweep. The edge sharpness cis originates from 
the beam pointing error of < 0.1 aiat and from variations 
in the magnetic bias field. The latter causes frequency 
fluctuations of ~ 5 kHz, which translate into an effective 
pointing error of 0.05 aiat at the maximum slope of the 
addressing beam profile. The resolution a a could in prin- 
ciple be further reduced by a narrower microwave sweep, 
at the cost of a larger sensitivity to the magnetic field 
fluctuations. A larger addressing beam power would re- 
duce this sensitivity, but we observed that this deformed 
the lattice potential due to the imperfect a" -polarization, 
allowing neighbouring atoms to tunnel to the addressed 
sites. 



Coherent tunneling dynamics 

The preparation of an arbitrary atom distribution 
opens up new possibilities for exploring coherent quan- 
tum dynamics at the single-atom level. As an example, 
we studied the tunnelling dynamics in a one-dimensional 
lattice (see Fig.[4| which allowed us to determine how 
much our addressing scheme affects the vibrational state 
of the atoms. We started by preparing a single line of 
up to 18 atoms along the y direction before we lowered 
the lattice along the x direction to = 5.0{5)Er within 
200 /is. At the same time, the other lattices were lowered 
to Vy = 30^^^ and Vz = 23 Er, which reduced the exter- 
nal confinement along the x direction, but still suppressed 
tunnelling in the y and z directions. After a varying hold 
time t, allowing the atoms to tunnel along x, the atomic 
distribution was frozen by a rapid 100 /is ramp of all lat- 
tice axes to 70 Er. By averaging the resulting atomic 
distribution along the y direction and repeating the ex- 
periment several times, we obtained the probability dis- 
tribution of finding an atom at the different lattice sites 
(Fig.|4j bottom row). 

This probability distribution samples the single-atom 
wave function after a coherent tunnelling evolution. We 
observed how the wave function expands in the lattice 
and how the interference of different paths leads to dis- 
tinct maxima and minima in the distribution, leaving 
for example almost no atoms at the initial position after 
a single tunnelling time (Fig.jl]^). This behaviour differs 
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FIG. 4: Ground state tunnelling dynamics, a, Atoms were prepared in a single line along the y direction before the lattice 
along the x axis was lowered, allowing the atoms to tunnel in this direction (b-d). The top row shows snapshots of the atomic 
distribution after different hold times. White circles indicate the lattice sites at which the atoms were prepared (not all sites 
initially contained an atom). The bottom row shows the respective position distribution obtained from an average over 10 — 20 
of such pictures, the error bars give the la statistical uncertainty. A single fit to all distributions recorded at different hold 
times (red curve) yields a tunnelling coupling of J^^^ /h — 940(20) Hz, a trap frequency of C(;trap/(27r) = 103(4) Hz and a trap 
offset of XofFs = -6.3(6) aiat- 



markedly from the evolution in free space, where a Gaus- 
sian wave packet disperses without changing its shape, 
always preserving a maximum probability in the center. 
For longer hold times, an asymmetry in the spatial dis- 
tribution becomes apparent (Fig.|4]i), which originates 
from an offset between the bottom of the external har- 
monic confinement and the initial position of the atoms. 
We describe the observed tunnelling dynamics by a sim- 
ple Hamiltonian including the tunnel coupling J^^^ be- 
tween two neighbouring sites and an external harmonic 
confinement, parameterized by the trap frequency cJtrap, 
and the position offset Xoffs (see Appendix). A single fit 
to all probability distributions recorded at different hold 
times yields J^^^ /h = 940(20) Hz, cjtrap/(27r) = 103(4) Hz 
and Xoffs = — 6.3(6) aiat- This is in agreement with the 
trap frequency a;trap/(27r) = 107(2) Hz obtained from an 
independent measurement via excitation of the dipole 
mode without the x lattice, whose contribution to the 
external confinement is negligible compared to the other 
two axes. From J^^\ we calculated a lattice depth of 
Vx = 4.6(1) which agrees with an independent cal- 
ibration via parametric heating. The expansion of the 
wave packet can also be understood by writing the ini- 
tial localized wave function as a superposition of all 
Bloch waves of quasi-momentum with — 7r/aiat < q < 
Tr/aiat- To each quasi-momentum hq, one can assign a 
velocity Vq = determined by the dispersion rela- 

tion E{q) = — 2J^^^ cos(gaiat) of the lowest band. The 
edges of the wave packet propagate with the largest oc- 
curring velocity I'max ~ 2J^°^aiatA = 1.88(4) aiat/ms, in 



agreement with our data. 

Our measurements constitute the first observation of 
the ground state tunnelling dynamics of massive parti- 
cles on a lattice with single-site resolution. A similar 
continuous-time quantum walk has been demonstrated 
with photons in an array of evanescently coupled pho- 
tonic waveguides [37]. For massive particles, a discrete 
quantum walk of single atoms has been observed us- 
ing a sequence of spin manipulations and spin-dependent 
transports in an optical lattice [38] and also with trapped 
ions [39 . Without single-particle and single-site resolu- 
tion, a continuous-time quantum walk in the ground state 
has been observed for ultracold fermionic atoms by mea- 
suring their ballistic expansion in a weak lattice |40J . 

In a second tunnelling experiment, we observed the 
faster dynamics of atoms in the first excited band (see 
Fig.[5|. For this, we deliberately excited the atoms by 
introducing a pointing offset 5x of the addressing beam, 
which caused a shift of the potential wells during the 
switch-on. We repeated the same tunnelling experiment 
as above with a hold time of t = 1 ms for different point- 
ing offsets 6x. For a small pointing offset {Sx = O.laiat 
in Fig.jsjD) we observed a narrow distribution, compared 
to a much broader one for a large offset {Sx = — 0.4aiat 
in Fig.[5|i). We attribute this to different fractions / of 
atoms in the first band which is characterized by the 
higher tunnelling rate J*^^^. We fitted the distribution of 
Fig.jS^i to a two-band model (see Appendix) and found 
J(i)/n = 6.22(6) kHz. This is in excellent agreement with 
the expected value of J^^^ /h = 6.14(6) kHz from a band 
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FIG. 5: Tunnelling dynamics of the first excited band. 

Some atoms were excited to the first band by a pointing offset 
5x of the addressing beam, a and b show the atomic position 
distribution after 1 ms tunneUing time for 5x = — 0.4aiat and 
5x — O.laiat, respectively. We fitted the data with a model 
that includes atoms in the zeroth band (red line) and a frac- 
tion / in the first band (orange line) (see Appendix). The 
right insets of a, b show corresponding original images. The 
left inset of b shows / versus 5x with a broad minimum of 
/ = 13(2)%, indicating that most of the atoms are left in the 
ground state. 

structure calculation in which we used J^^^ as an input 
parameter to calculate the lattice depth. Our measure- 
ment of the fraction of excited atoms / as a function of 
the pointing offset 5x (inset in Fig.jsjD) shows that the 
atoms are strongly heated for large pointing offsets. By 
contrast, only a small fraction of the atoms is excited to 
the first band for small pointing offsets \5x\ < O.laiat, 
yielding a ground state population of 1 — / = 87(2)%. 



Discussion 

In summary, we have demonstrated full two- 
dimensional single-site and single-atom spin control in 
an optical lattice with sub-diffraction-limited spatial res- 
olution. Starting from a Mott insulator with unity filling, 
we achieved a spin- flip fidelity of 95(2)% with negligible 
influence on the neighbouring lattice sites. Our scheme 
leaves most of the atoms in the vibrational ground state. 
The control of single spins in a strongly correlated many- 



body system on a lattice opens many new possibilities for 
studying quantum dynamics and quantum phases. Our 
technique will allow us to create out-of-equilibrium states 
or local perturbations in order to observe the ensuing 
dynamics of the many-body system, such as spin-charge 
separation [H |6] or spin impurity dynamics beyond Lut- 
tinger liquid theory [4 . Our studies of the tunnelling 
dynamics at the single-atom level can be extended to cor- 
related particle tunnelling j4lT[43] , also in higher dimen- 
sions, or to observe transport across local impurities [44] 
or potential barriers. The tunneling can also be used to 
circumvent the pairwise losses during the imaging [27l [28] 
by letting the atoms of a one-dimensional system spread 
along the perpendicular direction in order to obtain a 
sufficiently low density. Further prospects are the imple- 
mentation of novel cooling schemes relying on the local 
removal of regions with high entropy [8 , 45 . The single- 
spin control in our large systems with several hundreds 
of atoms also opens new perspectives for scalable quan- 
tum computing. Combining single-qubit manipulation 
with local readout and a global entanglement operation 
in a spin-dependent lattice ^461 147j would be the basis of 
a one-way quantum computer [TTJ [12] . For the circuit 
model of a quantum computer, two-qubit operations can 
be realised by Rydberg gates between selected atom pairs 
in the lattice [48, |49^ . 
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Appendix 



Calibration of the addressing beam position 

To move the addressing laser beam in the object plane, 
we changed the angle of the beam entering from the re- 
verse direction into the microscope objective using a two- 
axis piezo mirror. The device has an angular resolution 
of 5 /irad, corresponding to a theoretical resolution in the 
object plane of 0.02 aiat — 10 nm. In order to position 
the addressing laser beam onto the atoms with high pre- 
cision, we measured calibration functions that translate 
the two control voltages of the piezo mirror into image 
coordinates. This calibration was performed by replacing 
the far detuned addressing laser beam by a near resonant 
molasses beam that follows the identical beam path. Us- 
ing in addition the x and y molasses beams, we took a 
fluorescence image of a large thermal atom cloud in the 
vertical lattice alone and observed a strongly enhanced 
signal at the position of the focused beam. We deter- 
mined the position of this fluorescence maximum with 
an uncertainty of 0.2 pixels in our images, corresponding 
to 0.05 aiat = 25 nm in the object plane. The long term 
drifts of the addressing beam position are on the order of 
0.1 aiat per hour, which we took into account by regular 
recalibration of the beam position. 
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Lattice phase feedback 

In order to compensate slow phase drifts of the opti- 
cal lattice, we applied a feedback on the position of the 
addressing beam. We determined the two lattice phases 
along X and y after each realisation of the experiment 
by fitting the position of isolated atoms. Averaging over 
the positions of typically 1 — 5 isolated atoms per image 
allowed us to determine the lattice phase to better than 
0.05 aiat- For the determination of the phase, we used the 
lattice constant and the lattice angles determined from a 
fluorescence image with many isolated atoms [28]. Since 
our phase drifts were slower than 0.04 aiat between two 
successive realisations of the experiment (25 s cycle time), 
we used the lattice phase from the last image to correct 
the addressing beam position. This feedback was done by 
adding appropriate offsets to the piezo control voltages. 



Microwave sweeps 

Our microwave sweeps are HSl-pulses [50] with time- 
dependent Rabi frequency Q{t) and detuning S{t) given 

by 



^}{t) = sech 
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d{t) 



tanh 



2t 



(la) 



(lb) 



where Qo/{27t) = 3 kHz is the maximum Rabi-frequency, 
/3 = 5 is a truncation factor, Tp = 20 ms is the pulse 
length, and crMw/(27r) = 60 kHz is the sweep width. The 
detuning d{t) is measured relative to the center of the 
sweep at cjmw = — Amw (see Fig.jljD). Here, cjq is the 
bare resonance between the two hyperflne states, includ- 
ing the shift of -570 kHz due to the magnetic bias fleld 
along the z direction and AMw/(27r) = —75 kHz denotes 
the offset of the sweep center. 



Spin-flip fldelity 

In order to determine the spin-flip fldelity, we fltted the 
hole probability po as a function of the pointing offset Sx 
(see Fig.lsb) to a flat-top model function given by 



PoiSx) = I [erf (^^^) + erf + B. (2) 

Here, erf(x) = 2/y^ J^^ e^'^^dr is the error function, 
aa denotes the full-width at half-maximum of the flat- 
top proflle and cTg the edge sharpness. We chose this 
model function since our HSl-pulses (see above) produce 
a flat-top population transfer proflle, the edges of which 
are dominated by randomly fluctuating quantities (beam 
pointing and magnetic flelds) following Gaussian statis- 
tics. The addressing fldelity is deflned diS F = A/{1 — B) 
taking into account that the maximum hole probability 



These yield a hole with probability B at unsuccessfully 
addressed sites which occur with probability 1 — F, such 
that pg^^^ = F + (1 - F)B. 



Single-particle tunneling dynamics 

We describe the coherent tunneling dynamics on k 
2n + 1 lattice sites by the Hamiltonian 



= -J(0) ^ at a,_,) 

i = —n 
n 

+K3xt ^ {i-XoEsfafai, (3) 

i = —n 

where J^^^ is the tunnel coupling in the lowest band, 
af {di) is the creation (annihilation) operator for a 
particle at site i. The strength of the external har- 
monic potential with trapping frequency cjtrap is given 
by Kxt = ^^^trap<^Lt' ^offs describes a position off- 
set with respect to the bottom of the harmonic potential. 
The single particle wave function and its coherent time 
evolution are given by 



i = —n 

= exp 



(4) 



with the initial condition ^^^^(0) = ag" |0) and the vac- 
uum state |0). The resulting probability of flnding the 

particle at lattice site i after time t is P^^\t) = |c-^^(t)p. 
For analyzing the data of Fig. [4] we calculated the time 
evolution for k = 17 lattice sites. 



Tunneling in the first excited band 

When including tunnelling in the flrst band, we assume 
an incoherent sum P^°\t) of the distributions p(°)(t) of 
the zeroth and P^^\t) of the first band as 



pt°t(t) = (1 - /)p(0)(t) + fP^^\t). 



(5) 



Po 



A -\- B also includes holes from thermal defects. 



The Hamiltonian H^^"^ in the flrst band and the coher- 
ent dynamics are identical to the ones of the zeroth band 
(Eqs.||and|4|, except for a different tunnel coupling J*^^^. 
When fltting this model to our data, we kept entrap, ^offs 
and J*^^^ flxed at the values obtained from the data dis- 
played in Fig. [4] We extracted J^^^ from the data of 
Fig. [5^ and used this value to flt the results for other 
pointing offsets. For the data in Fig.|5) the parameters 
of our microwave sweep were such that also neighbour- 
ing atoms were addressed. We took this into account by 
summing over two distinct probability distributions with 
a second starting position in the direction of the pointing 
offset. 



